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Abstract 

Motivated by the very recent measurements performed at the LHCb and the Tevatron of 
the Bg — Bg mixing, in this paper we revisit it in a family non-universal Z' model, to 
check if a simultaneous explanation for all the mixing observables, especially for the like- 
sign dimuon charge asymmetry observed by the DO collaboration, could be made in such a 
specific model. In the first scenario where the Z' boson contributes only to the off-diagonal 
element M( 2 , it is found that, once the combined constraints from AM S , (p s and Ar s are 
imposed, the model could not explain the measured flavour-specific CP asymmetry a s j s , 
at least within its la ranges. In the second scenario where the NP contributes also to the 
absorptive part T s l2 via tree-level Z'-induced b — > ccs operators, we find that, with the 
constraints from AM S , 4> s and the indirect CP asymmetry in Bj — > J/^Ks taken into 
account, the present measured la experimental ranges for a s j s could not be reproduced 
too. Thus, such a specific Z' model with our specific assumptions could not simultaneously 
reconcile all the present data on B® — B® mixing. Future improved measurements from 
the LHCb and the proposed superB experiments, especially of the flavour-specific CP 
asymmetries, are expected to shed light on the issue. 
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1 Introduction 



Within the Standard Model (SM), the flavour- changing neutral current (FCNC) processes occur 
only at the loop level, and are therefore a very sensitive probe of new physics (NP) beyond the 
SM. In this direction, an outstanding role is played by the B® — B° s mixing, the phenomenon of 
which is described by two off-diagonal elements, M( 2 of the mass and rf 2 of the decay matrix. 
These two complex parameters can be fully determined by the following four observables: the 
mass difference AM S , the decay width difference Ar s , the CP- violating phase <f) s , as well as the 
flavour-specific CP asymmetry a s j s in £? s -meson decays [Tj [2] . 

The mass difference AM S was first measured at the Tevatron [3j H] , with the most precise 
published value given by the CDF collaboration [4] 

AM S CDF = (17.77 ± 0.10(stat.) ± 0.07(syst.)) ps^ 1 . (1) 

Employing a data sample of 340pb _1 taken in 2011, the LHCb collaboration has recently 
updated its analysis made in 2010 [5], and the most precise value of AM S is found to be [6] 

AM LHCb = ( 17 725 ± o.041(stat.) ± 0.026(syst.)) ps" 1 . (2) 

Both of these measurements agree quite well with the corresponding SM prediction El E] 

AMf = (17.3 ± 2.6) ps" 1 . (3) 

The decay width difference Ar s and the CP-violating phase S can be simultaneously deter- 
mined from an angular analysis of the decay B s — > J/ip<p. Both the CDF [TU] and DO [UJ col- 
laborations have recently updated their previous analysis based on 1.35fb _1 [12J and 2.8fb~ 1 [13J 
of data sample, respectively. While the two measurements agree within uncertainties, the CDF 
result constrains <p s to a narrower region, and we therefore quote [TO] 

0^ DF G [-1.04,-0.04] U [-3.10,-2.16] rad, 
Ar^ DF = (0.075 ± 0.035(stat.) ± 0.006(syst.)) ps" 1 , (4) 

given at 68% confidence level (C.L.). It is noted that these measurements have been superseded 
by the recent LHCb analysis [14] . Combining the two channels B s — > J/ip<f) [14] and B s — > 
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J/i^fo [IS], the LHCb collaboration presents the most precise measurement at 68% C.L. [T5] 

0LHCb = ^ 0? ± .i7(stat.) ± 0.06(syst.)) rad , 
Ar LHCb = 123 ± o.029(stat.) ± 0.011(syst.)) ps" 1 . (5) 
Compared to the corresponding SM predictions [7J [8j [9] 

^ M = -2&xg[-{V ts V? b )/{V cs VX)\ = (-0.037 ±0.002) rad, 
Arf* = (0.087 ± 0.021) ps" 1 , (6) 

these improved measurements show a better agreement with the SM expectations; especially, 
no evidence for a NP phase in — mixing is found by the LHCb collaboration [HJ H5] . 

The flavour-specific CP asymmetry a s j s can be extracted from a measurement of the like-sign 
dimuon charge asymmetry A b sl of semileptonic b-hadron decays [IB]. Employing a data sample 
of 9.0ffcT 1 of pp collisions, the DO collaboration has recently updated its previous analysis [17] . 
and the new measurement reads [T8] 

a j,do = (_ 7g7 ± .172(stat.) ± 0.093(syst.))% , 

a s f f° = (-1.81 ±1.06)%, (7) 
which, compared to the corresponding SM predictions [7J El |9] 

A b f M = (-2.0 ±0.3) x 10~ 4 , 

a s,SM = (Lg ± Q 3) x 10 - 5 _ (g) 

differ by about 3.9<r and 1.7cr, respectively. 

Motivated by the above observations, it is interesting to investigate if a specific NP model, 
while satisfy the constraints from AM S and <p s , could also simultaneously explain the measured 
values of Ar s and a s j s . It is also expected that the combined constraints from these observables 
could provide further information about the model parameter space. 

As is well-known, a Z' boson with family non- universal, flavour-changing couplings could 
arise in many well-motivated extensions of the SM with an additional U(l)' gauge symmetry [TjJJ 
|2"U| |2"T] . Searching for such an extra Z' boson is an important mission of the Tevatron [22j and 
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LHC [23] experiments. Performing constraints on the new Z' couplings through low-energy 
precise processes is, on the other hand, very crucial and complementary for these direct searches. 
It is interesting to note that such a family non-universal Z' model could bring new CP-violating 
phases beyond the SM and have a large effect on many FCNC processes [HI [20j [21] , such as the 
B s — B s mixing [2H [2S1 I2E1 123 I2E] j as well as some rare [2H] and hadronic B-meson decays [30J . 

Thus, in this paper we shall revisit the B s — B s mixing in a family non-universal Z' model, 
and check if it could simultaneously explain the measured values of AM S , AT S , S , as well as a^ s , 
updating our previous analysis made in Ref. [21], where only effects of the flavour-changing left- 
handed Z' couplings, with all the right-handed couplings being flavour-diagonal, on AM S and 4> s 
were considered. Furthermore, motivated by the like-sign dimuon charge asymmetry observed 
by the DO collaboration [18], we shall also consider a scenario with sizable Z' contribution to 
the off-diagonal element r\ 2 coming from the four-quark operators of the form {sb)v - a{cc)v±a- 

Our paper is organized as follows. In Sec. [2j we first recapitulate the theoretical framework 
for B® — B® mixing, and then discuss the Z'-boson contribution. In Sec. [3} we give our detailed 
numerical results and discussions. Our conclusions are made in Sec. [4] The relevant input 
parameters are collected in the appendix. 



2 Theoretical framework for B® — B® mixing 

In this section, we shall first recapitulate the theoretical framework for B® — B° s mixing within 
the SM. The family non-universal Z' model, as well as its contribution to the off-diagonal 
elements M( 2 and r^ 2 , which govern all the observables in B® — B® mixing, are then discussed. 



2.1 Observables in B). — B^ mixing 

The phenomenon of B® — B® mixing is described by a Schrodinger equation [Tj [9] 

with the mass matrix M s = M s ' and the decay matrix V s = . By diagonalizing M s — |r s , 
one can obtain the two mass eigenstates \B L ), \B H ) with masses M L , M H and decay rates T L , 
Th- Here L and H indicate the light and the heavy state, respectively. The mass and the width 
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difference between the two states \B H ) and \B L ) are then defined, respectively, as [H |9] 

AM S = M H -M L = 2\M* 2 \ , 

Ar s = r^-Tir = 2|r; 2 |cos$ a , (10) 

where $ s = arg(— Mf 2 /rf 2 ) is the CP-violating phase, and numerically irrelevant corrections 
that are proportional to (rf 2 /M;f 2 ) 2 have been neglected here. 

The CP-violating phase cj) s appears in b — > ccs decays of B s meson, taking possible mixing 
effect into account. Taking the decay B s — > J/ip<p as an example, since it is dominated by the 
tree-level b — > ccs transition that is real within the SM, the measured CP-violating phase (f) s is 
a direct probe of the phase of B® — B° s mixing. Assuming that there are no NP contributions 
to the decay amplitude^, we can therefore write jH] 

s = argMf 2 = ^M + NP ) (n) 

with the second term denoting the NP phase contributing to the off-diagonal element M[ 2 . 
In terms of the two off-diagonal elements M( 2 and r^ 2 , the flavour-specific CP asymmetry 



a s j s can be written as 



V s IF 
1 12 \ _ I 1 121 



a s f = Im(-^) = l~j±\sm$ s . (12) 

As the semileptonic decay B s — > X£i>£ is a typical example of flavour-specific decays, this 
asymmetry is usually also named as the semileptonic CP asymmetry [2]. 

Thus, in order to predict the mixing observables AM S , AT S , <p s , as well as aj s , we need to 
know the off-diagonal elements M( 2 and T\ 2 both within the SM and in the Z' model. 

2.2 The family non-universal Z' model 

In many well-motivated NP models with an additional U(l) gauge symmetry, such as string 
constructions and/or grand unified theories, the associated Z' boson can generally have family 
non-universal, flavour-changing couplings [19]. The general formalism of such a model has been 
detailed in Refs. 



X NP effects that contribute to the off-diagonal element Tf 2 might also give a contribution to the decay [27l 
I [3T] . For convenience, in this paper we shall not consider this case for B s meson. 



In the physical mass-eigenstate basis, due to the non-diagonal chiral coupling matrix, the 
tree-level FCNC interactions appear both in the left- and in the right-handed sectors [20, l2~Tj . 
The relevant interaction Lagrangian can be written as [28] 

C Z > = -\- [{B^ b rfP L b + B« s^P R b + h.c.) + {B L qq q^P L q + B R q^Pr q)] % , (13) 

COS t/p^/ 

where Pl,r = |(1 =F 75) project onto the left- and right-handed chiral fields, g is the SM 
SU(2) coupling, and 9w the weak mixing angle. For convenience, we have absorbed the U(l)' 
coupling constant into the factors B { - , and written the above Lagrangian in terms of the SM 
coupling. It should be noted that, while the flavour- changing couplings are in general complex, 
the diagonal ones have to be real due to the hermiticity of the Lagrangian. 



Starting with the Lagrangian Eq. (13), and integrating out the heavy Z' boson, one can 
easily obtain the resulting effective |A£?| = 1 and |A£?| = 2 four-fermion interactions induced 
by tree-level Z' exchange 



2.3 The off-diagonal element M( 2 

The off-diagonal element M( 2 can be decomposed into the SM and Z' parts 

M{ 2 = (M* 2 ) SM + (M; 2 ) z > , (14) 

where the SM contribution could be found, for example, in Refs. [331 In our case, due to 
the simultaneous presence of left- and right-handed currents, the complete set of dimension-six 
operators consists of the following ones [36], |37] 

qVll = (s a lfM P L b a )(s?rPLbP) , QX RR = (s a lfM P R b a )(s ls YPRb' S ) , 

Q^ R = (s a ^P L b a )(s^P R bP) , Q L 2 R = {rP L b a ){s?P R bP) , (15) 

where a and (3 denote the colour indices. Because of the V— A structure of W /± -boson exchanges, 
only Q\ LL contributes within the SM. The operators QY and Qf R are generated at the scale 
fiz> ~ tuz', whereas Q 2 R appear through renormalization group (RG) running to some scale 
different from \iz<- 

Normalized to the effective Hamiltonian for \AB\ = 2 transition within the SM [35] 

U™(\ AB\ = 2) = ^ ml (V tb V t * s y C^) Q\ LL {n) + h.c. , (16) 
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the NP effective Hamiltonian induced by tree-level Z ! exchange can be written as [2TI 137] 



Uf s (\AB\ = 2) 



r<2 



m* w (v tb v t :r Acr L (») <?r» + a<7™m q™» 



\VRR/ 



+ h.c. 



(17) 



where Gp is the Fermi coupling constant, and my/ denotes the VF-boson mass. The Wilson 
coefficients at the high scale \iz' are obtained by integrating out the heavy Z' boson, with the 
final results given explicitly as 

AC^ LL (fi Z/ 

ACf*(/^ 
AC 2 M (/^ 





64vr 2 


B sb B sb 


m\, 


\[2G F m 2 w 




m| 


64tt 2 


B R b B R 


m\, 


\/2G F m 2 w 




m| 


64tt 2 


IB sb B sb 






(VtbV t * s ) 2 







(18) 



where is the Z-boson mass. It is noted that the last coefficient ACj^ vanishes at the high 
scale /j,z i in the absence of QCD effects assumed by us. 

In order to include the QCD RG evolution during the calculation of M^ 2 , following the 
method proposed in Refs. [361 38], we shall evaluate the Wilson coefficients and the hadronic 
matrix elements of local operators at the high scale \iz>- The final Z'-boson contribution to the 
element M[ 2 can be written as 
1 



(M 



12) Z' 



2m b 
G 2 



(B° s \n z cS (AB = 2)\B° S ) 



48tt 



rn 2 w (V tb V*) 2 m B J 2 (ACf"(M + AC™W)) P^fa* 



LRi 



(19) 



where mg j and fs s are the mass and decay constant of B s meson, respectively. The coefficients 
Pi(n) collect compactly all RG effects from scales below \iz> and hadronic matrix elements 
obtained by lattice methods at low scales, and are defined as [36] 
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(20) 



Analytic formulae for -Pj Q (/i) given explicitly in terms of the QCD RG factors and the non- 
perturbative bag parameters could be found in Ref . [36J . 



2.4 The off-diagonal element Tf 2 

The off-diagonal element r* 2 , being an inclusive quantity stemming from decays into final states 
common to and B° s mesons, could be computed using the heavy quark expansion (HQE) [39, 
HQ], which is a simultaneous expansion in Aqcd/^6 and a; s (m&). It is now known to next- 
to-leading-order (NLO) both in AQ CD /m b J3U E2] and in a s (mb) [331 EH S3], with a recent 
theoretical update made in Ref. [32] . Motivated by the like-sign dimuon charge asymmetry 
observed by the DO collaboration [XT] [18]. the possibility to have sizable NP effects in rf 2 has 
been discussed in the literature [261 EH [28], HI] . 

Following the notation specified by Beneke et al. [331 HI], the off-diagonal element r* 2 , which 
is related via optical theorem to the absorptive part of the forward-scattering amplitude, can 
be written as 

with the transition operator T defined by 

T = \mi J d 4 xT[n c s(x)n cH (0)}. (22) 

Here "H e g- is the low-energy effective Hamiltonian for \AB\ = 1 transition, which is obtained 
by integrating out the heavy particles above the scale /if, ~ m b , such as the top quark, W 
and Z bosons within the SM, as well as the heavy Z' boson in our case. It can be generally 
decomposed as 

H eS = n™(\AB\ = 1) + H! S (\AB\ = 1) , (23) 

where a detailed review and explicit expressions of the SM part H^(\AB\ = 1) could be found, 
for example, in Ref. [35] . 

The Z' contribution to T s 12 arises from dimension-six operators of the form b — > sff (with 
/ denoting a light fermion), most of which are, however, constrained to be small by the rare 
B-meson decays. For example, the semileptonic operators b — > se + e~ , b — > sfi + fx~ are severely 
constrained by the b — > s£ + £~ processes [2H] , while the light-quark operators b — > suu, b —> sdd, 
b — > sss are strongly bounded by the various hadronic B-meson decays like B — > ttK and 
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B — > <f>K [30]. Furthermore, it is found that the contribution from b — > sr + r~ clearly fail 
to describe the B s — B s mixing data within 68% C.L. [27]. Thus, in this paper we shall only 
consider the Z' contribution coming from the four-quark operator b — > sec. 



Starting from the general couplings of Eq. (13), we get four types of operators, (sb)v-A(cc)v±A 
and (sb)v+A(cc)v±Ai where the former two are already present in the SM. Since the chirality 
of the flavour-changing part for the latter two operators is flipped, to run their Wilson coef- 
ficients down to the lower scale fib ~ mb, we have to calculate the corresponding anomalous 
dimensional matrix. Furthermore, due to the presence of these new operator structures beyond 
the SM, the calculation of their contributions to T* 2 needs also to be extended. Since both of 
these calculations are non-trivial, for simplicity we shall restrict ourselves to the case where 
only purely left-handed flavour-changing B^ b coupling is present, which is usually assumed in 
the literature (2B1 [27J [2S]. For the flavour-diagonal part, on the other hand, both the left- and 
right-handed CKM-favored couplings B^ c and B^ c are taken into account. 

With these assumptions, the final Z'-boson contribution to the off-diagonal element T s l2 can 
be written as I2T1 



n! s (\AB\ = i) = K 6 ^( s -6)^ A (cc) v _ A + ^^( s -6) v ^(cc)^] 

v ^ m z' 

= -^= V tb V* [AC 3 Q 3 + AC 5 Q 5 ], (24) 

where (qq')v±A — 97^(1 =t 75)9') an d i n the second line we have normalized the result to the 
SM effective Hamiltonian [35J. Thus, with our assumptions, the Z' effects can be understood 
as corrections to the Wilson coefficients of SM QCD-penguin operators Q3 and Q5. The new 
Wilson coefficients at the high scale [iz> are given, respectively, as 

AO^zO = -77^ 4 , AC B M = ~^~ BfrZ • (25) 

V tb V ts m z , VtbVts m z' 

Since no new operator structures beyond the SM appear in our case, the remaining calculation 
of the element T\ 2 is the same as that within the SM, the details of which could be found in 
Refs. 



In the presence of the effective Hamiltonian Eq. (24), the indirect CP asymmetry in B^ — > 
J/ipKs is also affected due to the NP contribution to the decay b —> sec, and the effective 
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measured sin2/3 meas is given by [281 EI] 

sin2/3 meas = sin2/3 + 2|r|cos2/3sin0cos(5 a , (26) 

where <fi and 5 a are the NP weak and the SM strong phase, respectively. The parameter |r| 
denotes the magnitude of the ratio between NP and SM b — > sec decay amplitudes. In terms 
of the Z' couplings, it is given explicitly as 



9m 2 



2 \-\* n B sb( B cc + B cc 
m Z> V cbV cs 0,2 



(27) 



"Z' 

where a<i ~ 0.17 is the effective Wilson coefficient within the SM. Numerically we shall use the 
estimation |r| < 26.5% at la [2S] as a constraint on the model parameter space. 

For convenience, when performing the RG running between the scales /i^' and we shall 
neglect the small RG effects in the six-flavour theory throughout the paper. It should be noted 
that all the above discussions apply also to the B® — B® mixing, with the replacement s — > d. 



3 Numerical results and discussions 

Equipped with the formalism discussed in previous section, as well as the theoretical input 
parameters listed in the appendix, we proceed to present our numerical results and discussions 
in this section. 



3.1 The SM predictions and the experimental data 

Within the SM, our predictions for the B® — B® mixing observables are listed in Table [TJ The 
theoretical uncertainties are obtained by varying each input parameter within its respective 
range specified in the appendix and adding the individual uncertainty in quadrature. The 
corresponding experimental data, which has been summarized in the introduction section, is 
collected in the second column. For convenience, we also give in Table [T] some relevant mixing 
observables in the S^-meson system that will be used later. 

As can be seen from Table [TJ our result for the mass difference AM S agrees quite well the 
experimental data, which is therefore expected to give a strong constraint on the NP parameter 
space. For the CP-violating phase <p s , while the SM prediction is quite precise, the experimental 
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Table 1: Our theoretical predictions for the — B® mixing observables within the SM, with the 
corresponding experimental data given in the second column. For convenience, the SM predictions for 
some relevant i?9 — B®, mixing observables are also given. 



observable 



experimental data SM prediction 



AM S [ps- 1 ] 17.77 ±0.12 @] 17.27i|£j 

17.73 ±0.05 [6] 

<p s [rad] G [-1.04, -0.04] U [-3.10, -2.16] [TO] -0.0421±°;°°^ 

0.07 ±0.18 H5] 

Ar s [ps" 1 ] 0.075 ± 0.036 [TO] 0.065±g;^| 

0.123 ±0.031 [15J 

a s fs [%} -1.81 ± 1.06 |IH] 0.0033t°;°o?o 

AM d [ps -1 ] 0.507 ±0.004 [2] 0.5341°;^ 



[degree] 42.8 ± 1.6 [2] 48.57 



1.44 
-6.42 



data still has a large uncertainty, with the central value deviating from the SM result. Our 
prediction for the width difference Ar s , while being consistent with the CDF measurement, is 
about la below the LHCb measurement. There is, however, a tension for the flavour-specific 
CP asymmetry a s ^ s between the SM prediction and the DO measurement, which has triggered 
a lot of investigations both within the SM and in various NP models [261 EH EEJ HI] . Thus, it 
is interesting to check if the parameter space allowed by AM S could simultaneously explain the 
remaining observables, especially Ar s and a s j s in a specific NP model. 

In the following, we shall revisit the — B° s mixing in a family non-universal Z' model. 
We shall use two different data sets: the first one (Dl) corresponds to the first row and the 
second (D2) the second row listed in Table [TJ while for a s j s we use the latest DO measurement 
for both sets. In addition, for both data sets, our analyses are further divided into the following 
two scenarios: (Mf 2 )z/ ^ 0, (r* 2 )z' = (named as SI) and (Mf 2 )z' ^ 0, (rf 2 )z' 7^ (named as 
S2). To obtain the allowed parameter space, we shall use the theoretical predictions with la, 
as well as the experimental data with 2a error bars as constraint^} 

2 For the CP-violating phase <j> a , on the other hand, in the data set Dl we still use the 68% C.L. region 
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3.2 Z' effects in scenario SI 

In this subsection, we consider the first scenario SI with (M( 2 )z' 7^ and (rf 2 )z' — 0, where 
all the Z' effects are encoded in the flavour-changing couplings B^ b = \B^ b \e^ b and Bf b = 
I Bfb I e% ^ h • I n terms of the four parameters and (pg b R , and fixing the Z'-boson mass at 

mz> = 400 GeV, the off-diagonal element M( 2 can be formally written as 



M{ 2 = (Mi)sM 



1 + 2.11 IfifjlO- 2 ]! 2 e *(^+ 1 -2i°) + 2.11 |Bf 6 [10- 2 ]| 2 e 



16.41 |54[10- 2 ]| |52[10- 2 ]| e ^(^/2+^/2+1.21°) 



21 12 _2»(0«+1.21 o ) 



(28) 



where the magnitudes of the couplings |5 s6 ' | are given in unit of 10 . From Eq. (28), one can 



see that the fourth term has only a marginal effect if \B^ b \ ^> \Bf b \ or \B^ b \ <C \Bf b \, but gives 
the dominant contribution if \B^ b \ ~ \Bf b \. It is also noted that the Z'-boson contribution with 
a single left-handed coupling (the second term) is the same as that with a single right-handed 
coupling (the third term). 

Thus, based on the above observations, we shall consider the following two limiting cases: in 
case I Bg b varies arbitrary and Bf b = 0, while in case II the simplification B^ b = Bf b is assumed. 

3.2.1 Case I with Bf b arbitrary and Bf b — 

In this case, the Z' effect is parameterized by the two parameters \B^ b \ and which could 
be severely constrained by the measured observables AM S , <p s , AT S and a s j s . This is clearly 
shown in Fig. [TJ where the left and the right panel correspond to the data sets Dl and D2, 
respectively. The green and the blue areas correspond to the case by choosing a Z' boson 
with mz' = 400 GeV and mz> = 1 TeV, respectively, while the purple area is allowed by the 
measured a s j- s with la error bars. The second solution for <p s , <p s G [—3.10,-2.16] rad [ID] , is 
denoted by the pink area. 

From Fig. [I] we make the following observations: 

• as shown in the first two plots, decreasing the Z'-boson mass reduces the allowed pa- 
rameter space, which is also true for the other observables. Thus, for simplicity, in the 
following discussions, we shall choose a fixed Z'-boson mass with m^' = 400 GeV. 
measured by the CDF collaboration 10J . 
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Figure 1: The <^-|-^ibl parameter space in case I, allowed by the B® — B® mixing observables. 
The left and the right panel correspond to the data sets Dl and D2, respectively. The green and the 
blue areas correspond to the case by choosing a Z' boson with mz> = 400 GeV and mz' = 1 TeV, 
respectively, while the purple area is allowed by the measured a s j s with la error bars. The second 
solution for <j> s , (f) s £ [—3.10, —2.16] rad [10], is denoted by the pink area. See text for details. 
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• these mixing observables have a different and complementary dependence on the Z' cou- 
pling Bg b , and their combined constraints are very strong, as shown in the last two plots, 
where the parameter space allowed by the individual observable is reduced significantly. 

• the pink area in the third plot that is allowed by the second solution for the CP-violating 
phase <p s , (p s G [—3.10,-2.16] rad JTU], is already excluded, once constraints from the 
measured AM S and Ar s are taken into account. 

• comparing the left and the right panel, the data set D2 gives a more stringent constraint 
on the magnitude \B^ b \. However, at present no constraint on the weak phase § L sb could 
be obtained for both data sets from the measured observables AM S , S and Ar s . 

• as shown in the last two plots, there is no overlap between the parameter space allowed by 
a s j s with la error bars, as well as the one allowed simultaneously by the three observables 
AM S , <p s and Ar s . Thus, within such a specific Z' model, it is difficult to explain the 
observed flavour-specific CP asymmetry a s j s by the DO collaboration [T8] . 

3.2.2 Case II with the simplification Bf b = Bf b 

In this case, due to the simultaneous presence of left- and right-handed currents, the Z'-boson 
contribution to M[ 2 is dominated through the generation of the LR operators that renormalize 
strongly under QCD (36J E7]. With the simplification B^ b = Bf b , contributions from the VLL 
and VRR operators are completely canceled by that from the LR operators. This can be clearly 
seen from the numerical result given in Eq. (|28]). 

Corresponding to case II, we show in Fig. [2] the 4>^ b -\B^ b \ parameter space allowed by the 
measured B° s — B® mixing observables. The other captions are the same as in Fig. H In this 
case, the following observations could be made: 

• due to the presence of the LR operators, the mixing observables show a quite different 
dependence on the NP weak phase 0f b than in case I. The allowed ranges for the magnitude 
\B^ b \ are significantly reduced, being about half of that obtained in case I. 

• as in case I, with constraints from the measured AM S , <p s and Ar s taken into account, 
we still could not reproduce the present measured la experimental ranges of a s j s by the 
DO collaboration [TS] . 
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Figure 2: The ^fk-l-Bfhl parameter space in case II, allowed by the measured B® — B® mixing 
observables. The other captions are the same as in Fig. [TJ 

• in both case I and case II, the mixing observables (f) s , AT S and aj s , which are all related 
to the phase of M[ 2 , show a strong dependence on the NP weak phase <j)^ b . 
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Thus, combining the above two limiting cases, it is concluded that the scenario SI could not 
explain the measured like-sign dimuon charge asymmetry by the DO collaboration [TS], once 
constraints from AM S , <ft s and Ar s are taken into account. This motivates us to consider the 
second scenario S2 with (Mf 2 )z' 7^ and (r* 2 )z' 7^ 0. 



3.3 Z' effects in scenario S2 

In the scenario S2, the Z' boson contributes both to Mf 2 and to T\ 2 . With our assumptions 



discussed in Sec. 2.4 the Z' effects are now parameterized by the four parameters \B\ 



B^ c and B^ c . In terms of them, the off-diagonal element rf 2 can be formally written as 

^,+i.77°) , n m n I r l h n- 2 l I R« A<fa+w°) 



r i2 - ( r i2)sM 



0.033 |^ 6 [10- 2 ]| B^e« +1 ' 77 > + 0.010 \B%[lQ- 2 ]\ B*e l 



b\i Ysbi 



(29) 



where the Z'-boson mass is fixed at = 400 GeV. One can see that, depending on the 
relative strength of the couplings B^ c and B^ c1 as well as the weak phase 0^, the Z' boson 
could contribute either constructively or destructively to the off-diagonal element T s l2 - Since 
the flavour-changing couplings \B^ b \ and <pg b could be constrained by the measured observables 
AM S and <p s , we shall focus on the flavour-diagonal couplings B^ c and B^, which could be 
bounded by the indirect CP asymmetry in Bj — > J/ipKs [28, 31 J . 

Thus, in this subsection, imposing the constraints from AM S , <f> 3 , as well as the indirect CP 
asymmetry in Bd — > J/ipKs [28l[3T], we shall investigate if the specific Z' model could explain 
the measured Ar s and a S f s . For simplicity, with respect to the flavour-diagonal couplings B^ c 
and B^ c , our analyses are further divided into the following three different cases: 

• Case I with B^ c arbitrary and B^ c = 0. 

• Case II with B^ c arbitrary and B^ c = 0. 

• Case III with the simplification B^ c = B^ c . 

With the above assumptions, the Z' effects can be totally parameterized by the three parameters 
\Bg b \, (ftsb an d B^ c (or B^.). With the constraints from AM S , <p s and the indirect CP asymmetry 
in Bd — > J/ipKs imposed, our predictions for Ar s and a s j s in such a Z' model are shown 
in Figs. |3j [4] and |5j corresponding respectively to the three different cases listed above. In 
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these figures, the horizontal and vertical lines correspond to the experimental data with la 
error bars, while the SM prediction is indicated as a black point. The left and the right panel 
correspond to the data sets Dl and D2, respectively. Here we have fixed the Z'-boson mass at 
m z , = 400 GeV. 

From the correlation plots shown in Figs. [3j [4] and [5j the following observations are made: 

• in all these three different cases, once constraints from AM S , <p s and the indirect CP 
asymmetry in Bj — > J/ipKs are imposed, we still could not explain the measured like- 
sign dimuon charge asymmetry by the DO collaboration [18], at least within its la ranges. 
This is due to the fact that the flavour-diagonal couplings B^ c and B^ c are already severely 
constrained by the indirect CP asymmetry in B^ — > J /ipKg- 

• for the data set Dl, the predicted values of Ar s in these three different cases are all 
compatible with the experimental data within la error bars. However, for the data set 
D2, the result obtained in Case II is below the present measured la experimental ranges 
by the LHCb collaboration, with the other two cases being consistent with the data. 

• the flavour-specific CP asymmetry aj s shows a strong correlation with the observables 
AM S and <f) s , and it is therefore expected that more stringent information about the 
former could be obtained from the latter two observables, all of which will be measured 
more precisely by the LHCb collaboration in the near future. 

Thus, even with a contribution from tree-level Z'-induced b — > ccs operators included, a 
large correction to rf 2 in such a specific Z' model is excluded. It is therefore impossible to 
give a simultaneous explanation for the latest measurements of B® — B° s mixing observables, 
especially for the like-sign dimuon charge asymmetry observed by the DO collaboration |18j . 

3.4 Z' effects on B Q d — mixing 

For completeness, in this subsection we discuss the Z' effects on the B® — B® mixing observables 
AMd and <pd, assuming that the NP contribution comes only from the flavour-changing Z' 
couplings B^ R . As can be seen from Table [TJ our SM predictions and the experimental data 
for these two observables agree quite well with each other. It is therefore expected that these 
measurements could exert strong constraints on the Z' parameter space. 
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Figure 3: The correlation plot between AT S , a s j s and <p s , AM S for Case I where only the couplings 
B^ b and B^ c are present. The horizontal and vertical lines correspond to the experimental data with 
lo" error bars, while the SM prediction is indicated as a black point. The left and the right panel 
correspond to the data sets Dl and D2, respectively. The Z'-boson mass is fixed at mz' = 400 GeV. 
See text for details. 
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Figure 4: The same as in Fig. [3] but for Case II where only the couplings B^ b and B^ are present. 
The other captions are the same as in Fig. [3] 



Following the same procedure as in — B° s mixing, our final results are shown in Fig. [6j 
where the left and the right panel correspond to the case with B% b arbitrary and B^ b = (the 
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Figure 5: The same as in Fig.[i]but for Case III where both the couplings B^ b and B^ R are present, 
with the simplification B^ c = B^ c . The other captions are the same as in Fig. [3] 
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Figure 6: The <i>db'\^db\ parameter space allowed by the B® 
The left and the right panel correspond to the case with B db arbitrary and B db = 0, as well as the 
case with the simplification B db = B^\, respectively. The green and the blue areas are obtained by 
choosing a Z' boson with mz' = 400 GeV and mz' = 1 TeV, respectively. 

a Z' boson with rriz> = 400 GeV and mz' = 1 TeV, respectively. 

As can be seen from Fig. [6j the allowed parameter space in the 4>db~\-^db\ plane is greatly 
reduced with the combined constraints from AMd and <ftd imposed. It is also interesting to note 
that, similar to the hierarchy of the CKM matrix elements \(V t bV t * d ) / (V t bV t *)\ ~ 0(10^), the 
final allowed magnitudes l-B^I and l-B^I satisfy roughly the same hierarchy, \B db ' \/\B sb ' | ~ 
O(l0~ l ), within such a family non-universal Z' model. This may imply that the flavour- 
changing Z' couplings are also linked to the known structure of the SM Yukawa couplings [24J . 
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4 Conclusions 



In this paper, motivated by the very recent measurements performed at the LHCb and the 
Tevatron of the mass difference AM S , the decay width difference Ar s , the CP- violating phase 
4>si as well as the like-sign dimuon charge asymmetry in semi-leptonic b-hadron decays, we 
have revisited the B° s — B° s mixing in a family non-universal Z' model, to check if the specific 
model could simultaneously explain the present measured values of these observables. Our 
main conclusions are summarized as follows: 

• In both data sets Dl and D2, the combination of measured mass difference AM S and 
CP-violating phase <p s could give a strong constraint on the flavour-changing Z 1 couplings 
B s b i especially in the case where both the left- and right-handed currents are present at 
the same time. 

• In the first scenario where the Z 1 boson contributes only to the off-diagonal element 
Mf 2 , we find that, once constraints from AM S) <fi s and Ar s are taken into account, the 
specific model could not reproduce the present measured value of the flavour-specific CP 
asymmetry a^ s within its la ranges. 

• Motivated by the failure of the first scenario, we have then considered the second scenario 
where the Z 1 boson contributes both to M{ 2 and to Y\ 2 . For simplicity, we have assumed 
that the NP contribution to rf 2 comes only from tree-level Z'-induced b — > ccs operators, 
with only the presence of left-handed flavour-changing coupling B^ b . 

• In the second scenario, due to the presence of flavour-diagonal couplings B^ R , we have 
considered the constraint from the indirect CP asymmetry in Bd — > J/ipKs, in addition to 
the ones from AM S and 4> s . With these constraints imposed, we find that a large correction 
to rf 2 is already excluded. While the predicted value of Ar s is roughly compatible with 
the experimental data, the measured lex experimental ranges of a s j- s still could not be 
reproduced in such a NP model. 

• For completeness, we have also presented the Z 1 effects on the mixing observables AM d 
and (pa, assuming that the NP contribution comes only from the flavour-changing cou- 
plings B^ R . It is interesting to note that, similar to the hierarchy of the CKM matrix 
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elements, |(V^V^)/(V^ 6 V;*)| ~ C^IO" 1 ), the allowed magnitudes \B% R \ and \B^ R \ sat- 
isfy roughly the same hierarchy, \B^ R \/\B^ R \ ~ (9(10 _1 ). This may imply that the 
flavour-changing Z' couplings are also linked to the known structure of the SM Yukawa 
couplings [21]. 

In conclusion, the specific Z' model we are considering could not simultaneously explain 
the measured B® — B® mixing observables, especially the like-sign dimuon charge asymmetry 
observed by the DO collaboration. Future improved measurements from the LHCb and the 
proposed superB experiments, especially of the flavour-specific CP asymmetries, are expected 
to shed light on the issue. 
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Appendix: Theoretical input parameters 

In this appendix, we collect all the relevant input parameters for the B® — B® mixing observables 
both within the SM and in the family non-universal Z ! model. 

• The basic SM parameters 

First, we need some basic SM parameters, which are, if not stated otherwise, taken from the 
Particle Data Group [15] 

a s (m z ) = 0.1184 ± 0.0007, G F = 1.16637 x 10~ 5 GeV" 2 , sin 2 6 W = 0.23146, 

m z = 91.1876 GeV, m w = 80.399 GeV, m t = 173.2 ± 0.9 GeV [ffi], 

m Bd = 5279.50 MeV, m Bs = 5366.3 MeV, (30) 

where m t is the top-quark pole mass, and we use two-loop running for a s throughout this paper. 
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• The CKM matrix elements 

To discuss possible NP effects, we should treat the SM contribution as a theoretical background, 
and calibrate the CKM matrix elements exclusively on SM tree-level observables, which are 
largely insensitive to NP contributions. These include the moduli of CKM matrix elements from 
super-allowed (3 decays, leptonic and semileptonic meson decays, as well as the CP-violating 
phase angle 7 from tree-dominated B-meson decays [15]. By fitting these constraints to the 
Wolfenstein parametrization [17] of CKM matrix up to C(A 4 ), it is found that [48J 

A = 0.799 ± 0.026 , A = 0.22538 ± 0.00065 , p = 0.124 ± 0.070 , r) = 0.407 ± 0.052 , (31) 

which are used through this paper. 

• The heavy- and light-quark masses 

For the light-quark masses, we adopt the values determined by the Flavianet Lattice Averaging 
Group (FLAG) [UJ 

m s (2 GeV) = 94 ± 3 MeV, m d (2 GeV) = 4.67 ± 0.20 MeV. (32) 

The scale-invariant b- and c-quark masses are taken as [15] 

m b (m b ) = 4.191°;^ GeV, m c (m c ) = 1.29±g;^ GeV. (33) 

To get the corresponding pole and running masses at different scales, we use the NLO MS-on- 
shell conversion and running formulae collected, for example, in Ref. [5D] . 

For the mass parameter m% ow , which appears in the 1 /mt,-suppressed matrix elements for 
the off-diagonal element rf 2 , we take m^ ow = 4.7 ± 0.1 GeV [32J. 

• The non-perturbative parameters 

For the B-meson decay constants, we use the averaged lattice results performed by the CKM- 
fitter group [5T] 

f Bs = 231 ± 3 ± 15 MeV, / B ,// fl<l = 1.235 ± 0.008 ± 0.033, (34) 

where the first error is statistical and accountable systematic, while the second stands for 
systematic theoretical uncertainties [5"T] . 
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Among the non-perturbative bag parameters Bf, the parameter B\ LL is already known 
from the SM analyses, and we use [9] 

B i,B,(Vb) = 0.841 ±0.013 ±0.020, = 1.01 ± 0.01 ± 0.03, (35) 

renormalized at the scale \i b = 4.6 GeV. Since the parameter BY RR is the same as bY ll , 
we can get the former from the latter. For the SM contribution, on the other hand, the RG 
invariant parameters B B are adopted, and numerically we use 



B Ba = 1.291 ± 0.025 ± 0.035, B B jB Bd = 1.024 ± 0.013 ± 0.015, (36) 

obtained by averaging over the lattice results. 

To obtain the coefficient Pi{fiz>), we need the bag parameters B^ R and B^, the values 
of which are taken as |52] 



£fj> 6 ) = 1.79 ±0.04 ±0.18, B L x %jB L x % d = 1.01 ±0.03, 

B^ s (fi b ) = 1.14 ± 0.03 ± 0.06, B l 2 %JB^% a = 1.01 ± 0.02. (37) 

To predict the off-diagonal element r^ 2 , we also need the following bag parameters P, [32] 

Bs(fi b ) = 0.91 ± 0.03 ± 0.12, B Ro i (/i 6 ) = B Ai2s (fx b ) = 1.0 ± 0.5, 

5 2 

B R2 (fi b ) = B k ^ b ), B Rs (fi b ) = -B R3 (n b ) + -B A2 (n b ). (38) 
• The mass of Z' boson 

Since we are mainly concerned with the Z' couplings to fermions, we shall specify the Z'-boson 
mass. At present, the best bound is provided by the CDF measurement |54j . which rules out a 
Z' boson with a mass below 399 GeV at 95% C.L.. As a guideline, we consider both a relatively 
light Z'-boson with mz> = 400 GeV, as well as a heavier one with mz> = 1 TeV. 
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